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Accurate  ground-state  calculations  of  H2''^  using  basis  sets  of  atom-centered  Slater-type  orbitals 

Herbert  W.  Jones  and  Babak  Etemadi 
Department  of  Physics,  Florida  A&M  University,  Tallahassee,  Florida  32307 
(Received  24  September  1992) 

The  variational  method  is  used  to  find  the  ground-state  energy  of  the  molecular  hydrogen  ion,  by 
employing  a  linear  combination  of  atomic  orbitals  (LCAO)  of  the  Slater  type,  all  with  the  same  exponen¬ 
tial  parameter  For  the  equilibrium  position,  R  =2  a.u.,  eight-digit  accuracy  was  achieved  using  21  or¬ 
bitals  on  each  nucleus,  with  the  highest  principal  quantum  number,  N=6.  Beyond  R  =50,  thirteen¬ 
digit  accuracy  is  obtained.  The  wave  function  is  given  for  various  LCAO  at  the  equilibrium  position. 

Strategies  developed  here  will  be  useful  in  our  studies  of  polyatomic  nonlinear  molecules. 

PACS  number(s):  31.20.Di 


INTRODUCTION 


LCAO  METHOD  FOR  Hj^ 


Basis-set  methods,  i.e.,  the  algebraic  approximation 
method,  in  molecular  electronic  structure  calculations 
have  been  shown  by  Wells  and  Wilson  [1]  to  give  excel¬ 
lent  results  for  simple  molecules,  includ’ng  They 
employed  an  even-tempered  basis  set  of  Gaussian-type 
orbitals  (GTO’s).  In  this  paper,  we  use  the  more  “intract¬ 
able”  [2]  Slater-type  orbitals  (STO’s)  as  basis  with  the 
hope  of  developing  strategies  so  that  we  will  be  able  to 
take  advantage  of  their  more  physical  nature  for  use  in 
scattering  phenomena  and  molecular  structure  calcula¬ 
tions.  Wells  and  Wilson  also  emphasize  the  advantages 
of  basis-set  methods  over  methods  employing  numerical 
grids,  especially  in  regard  to  nonlinear  polyatomic  mole¬ 
cules. 

Bishop  and  Schneider  [3],  in  their  study  of  transforms, 
remarked  that  completeness  in  basis  sets  for  s  orbitals 
may  be  achieved  by  variation  of  the  principal  quantum 
number  or  the  exponential  parameter.  Thinking  of  the 
first  possibility,  Jones,  Etemadi,  and  Brown  [4]  used  s  or¬ 
bitals  from  Is  to  8s  to  obtain  convergence  with  these 
spherically  symmetric  orbitals  on  the  nuclei  of  H2^.  En¬ 
couraged  by  these  results,  we  extend  our  study  by  use  of 
an  ordered  sequence  of  orbitals  with  increasing  principal 
quantum  numbers  and  angular  momentum.  (All  orbitals 
have  the  same  orbital  parameter  or  screening  constant  ^.) 
By  this  procedure,  the  ground-state  energy  of  H2''^  for 
nuclear  separations  from  0.2  to  SO  a.u.  is  calculated  to  an 
accuracy  of  from  six  to  thirteen  digits.  In  addition,  the 
wave  function  for  the  equilibrium  position  (R  =2.0  a.u.) 
is  given  for  various  numbers  of  orbitals. 

Finkelstein  and  Horowitz  [5]  were  the  first  to  use  is  or¬ 
bitals  on  each  nucleus  with  an  adjustable  orbital  ex¬ 
ponent.  Dickinson  [6]  (corrected  by  Geller  and  Frost  [7]) 
used  Is  and  2p  on  each  nucleus.  Miller  and  Lykos  [8] 
used  a  Is,  2p,  and  3d  orbital  on  each  nucleus.  (We  note 
that  our  result  for  a  Is,  2s,  and  2p  orbital  on  each  nucleus 
is  more  accurate.) 


In  the  variational  treatment  of  a  quantum-mechanical 
system,  the  rjt  function  with  one  or  more  adjustable  con¬ 
stants  or  parameters  is  assumed  as  the  solution  of  the 
problem,  and  then  the  parameters  are  adjusted  to  mini¬ 
mize  the  energy.  For  our  specific  example,  we  take  a 
linear  combination  of  atomic  orbitals  (LCAO),  all  with 
the  same  value  for  the  nonlinear  parameter  f,  for  our  tj) 
function.  A  STO  has  the  form 

X=Ar^-'e-^'Y^{e,(p),  A  = 

We  take  a  straightforward  linear  combination  of  STO’s, 

N-\ 

tl>=K  2  2  CsLA^r^-'e-^'Yl{e,,p)  , 

\  =  )  L=0 

as  our  LCAO  for  the  ground  state  (M  =  0)  and  center  a 
set  on  each  nucleus,  i.e.,  nucleus  a  and  b.  K  is  the  nor¬ 
malization  constant  chosen  such  that 

J’xp*Vl/dT=l  ('l'  =  tfr„+l/-fc)  . 


TABLE  I.  Exponential  parameter  f  at  the  lowest  energy  for 
various  values  of  the  intemuclear  distance  R  and  maximum 
principal  quantum  number  for  a  LCAO,  pertaining  to 

IV, _ 


^liigh 


R 

1 

2 

3 

4 

5 

6 

0.2 

1.94 

2.12 

2.50 

2.50 

1.0 

1.54 

1.91 

1.93 

2.11 

2.20 

2.30 

2.0 

1.24 

1.58 

1.70 

1.76 

1.93 

2.00 

3.0 

1.09 

1.34 

1.56 

1.57 

1.68 

1.80 

20.0 

1.00 

0.98 

0.98 

0.97 

1.00 

1.00 

50.0 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 
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^NL  are  the  linear  parameters  and  iVhigt,  is  the  value  of 
the  largest  principal  quantum  number  used  for  a  given 
trial  (it  will  vary  from  1  to  6). 

We  write  the  Schrodinger  equation 


—  Ov 
00  ^ 
as 
ro 

r4  fs 
O  'O 
00  m 
sc 
oo 

ri  m  o 
os  rf  ^ 


and  substitute  in  the  ^  under  consideration,  and  arrive  at 
a  secular  equation  for  the  energy  (eigenvalue)  [9].  This  ei¬ 
genvalue  problem  was  solved  on  a  CRAY-YMP  comput¬ 
er  using  an  IMSL  program.  The  matrix  elements  for  this 
problem  can  essentially  be  written  as  combinations  of 
overlap  integrals  [10].  Tai  [11]  pointed  out  numerous  er¬ 
rors  in  early  work  that  evaluated  overlap  integrals.  It  is 
only  recently  that  reliable  methods  for  the  evaluation  of 
STO  overlap  integrals  with  high  quantum  numbers  have 
been  developed  (see  Weniger  and  Steinbom  [12],  and 
Jones  [13]). 


RESULTS 

Table  I  shows  the  values  for  ^  that  give  the  minimum 
energy  for  various  values  of  the  intemuclear  distance  and 
^high-  Table  II  shows  the  corresponding  electronic  ener¬ 
gies,  The  exact  energy  for  the  first  row  was  taken  from 
Wind  [14].  The  other  exact  energies  were  taken  from 
Madsen  and  Peek  [IS].  (The  values  were  obtained  by 
subtracting  the  Coulomb  energy,  1/R,  of  the  nuclei  from 
the  electronic  energies.)  Using  A(high=6,  which  puts  21 
orbitals  on  each  nucleus,  and  R=2,  we  obtain  the  same 
energy  to  the  eight  digits  given  by  Wells  and  Wilson  [1] 
who  used  42  Gaussian-type  orbitals  on  each  nucleus.  The 
faster  convergence  of  STO’s  was  expected.  For  the  large 
separation  of  nuclei  at  /I  =  SO,  we  obtain  full  thirteen¬ 
digit  agreement  with  Madsen  and  Peek,  for  ^^high 
IS  orbitals  on  each  nuclei.  This  complete  agreement  is 
found  for  all  values  of  R  beyond  SO.  At  R  =0.2,  we  can 
only  use  ten  orbitals  ((Vhigh~'^J>  otherwise,  linear  depen¬ 
dencies  prevent  the  IMSL  program  from  functioning. 
Nevertheless,  even  at  the  small  distance  of  R  =0.2,  prac¬ 
tically  six-digit  accuracy  is  achieved.  Table  III  records 
the  orbital  coefficients  obtained  from  the  eigenvectors  of 
the  energy  minima  at  the  equilibrium  distance,  R  =2  a.u., 
for  a  various  number  of  orbitals  determined  by  N),ig|,. 


CONCLUSION 

We  have  shown  that  STO’s  are  a  viable  alternative  to 
OTO’s  when  dealing  with  a  simple  one-electron  molecule. 
Furthermore,  it  has  been  demonstrated  that  orbitals  with 
varying  principal  quantum  numbers  are  an  effective  re¬ 
placement  for  even-tempered  basis  sets  that  mostly  vary 
the  nonlinear  orbital  parameters.  The  Ldwdin  a-fuhction 
method  [16]  has  proven  to  be  effective  in  a  production 
mode  for  generating  accurate  values  of  overlap  integrals. 
More  complicated  molecules  are  under  investigation  us¬ 
ing  these  methods. 
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TABLE  III.  Coefficients  Cvi  before  orbitals  that  make  up  LCAO  with  highest  quantum  number 
^high.  at  the  equilibrium  distance  R  =2.0  a.u.,  for  H2^.  Each  nucleus  has  identical  orbitals.  K  is  the 
normalization  constant.  _  _ 


Orbital 

^ high  ~  2 

^high  ^ 

Arh.gh=4 

iVhi,h  =  5 

^high  ^ 

Is 

1.0 

1.0 

1.0 

1.0 

1.0 

2s 

0.597  502 

0.598  182 

0.636  512 

0.737  380 

0.745  909 

2pa 

0.185418 

0.130488 

0.112  303 

0.094421 

0.094957 

3s 

0.191573 

0.209  123 

0.448071 

0.400324 

3pa 

0.112262 

0.150186 

0.152  873 

0.153  366 

3da 

0.045904 

0.047706 

0.026401 

0.026  356 

4s 

0.060503 

0.076795 

0.164076 

4/>cr 

0.013  549 

0.069072 

0.076866 

4da 

0.013046 

0.051045 

0.053648 

Afa 

0.013  812 

0.016356 

0.016315 

5s 

0.111965 

0.031038 

ipa 

0.019  395 

0.023  338 

5d<T 

-0.005  287 

-0.000982 

5fa 

0.000791 

0.004847 

5ga 

0.003  948 

0.007663 

6s 

0.035  435  1 

6pa 

0.006673 

6d(T 

0.003  870 

6f<T 

-0.000462 

6ga 

-0.002  835 

6ha 

0.001  134 

K*~ 

0.131681 

0.105231 

0.093  3144 

0.063  955  1 

0.063  3508 
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